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The semi-relativistic approach to electron and neutrino quasielastic scattering from 
nuclei is extended to include final-state interactions. Starting with the usual non- 
relativistic continuum shell model, the problem is relativized by using the semi- 
relativistic expansion of the current in powers of the initial nucleon momentum and 
relativistic kinematics. Two different approaches are considered for the final-state 
interactions: the Smith- Wambach 2p-2h damping model and the Dirac-equation- 
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■ based potential extracted from a relativistic mean-field plus the Darwin factor. Using 

the latter the scaling properties of (e, e') and {y^pr) cross sections for intermediate 
momentum transfers are investigated. 

PACS numbers: 24.10.Jv; 21.60.Cs; 25.30.Fj; 25.30.Pt; 



I. INTRODUCTION 

In recent years, the so-called semi-relativistic (SR) approach has been explored as a 
convenient and easily implementable way to "relativize" existing Schrodinger-based models 
of quasielastic (QE) electron scattering from nuclei 

aaaana 

. This approach (which 



differs from other approaches also called "semi-relativistic" in the literature [3, S]) is based 
on the SR expansion of the on-shell electromagnetic current in powers of p/m^ only — i.e., 
the momentum of the initial nucleon divided by the nucleon mass — leaving untouched the 
dependence on the momentum transfer q and energy transfer u. Implementing this current 
expansion and employing relativistic kinematics, the resulting SR models can be applied 
at large values of the momentum transfer, of order a GeV or more, where the traditional 
relativistic corrections [9| based on expansions of the current in powers of l/ra^ are bound to 
fail, whatever the order of the expansion. The existence in the literature [lo| of expansions 
up to 4th order, at least in a Foldy-Wouthuysen transformation where they are seen to be 
quite complicated, raises simplicity as another advantage of the SR approach. In fact, by 
carefully grouping the expansion terms one can write them as (q, ^-dependent factors times 
the traditional [ll, [3] leading-order charge, convection and spin-magnetization current 
operators Moreover, the SR first-order (0(p/m^)) correction to the charge operator 



lv obtained in the 



131 ] is proportional to the spin-orbit charge operator iq ■ (cr x p), usua 
second-order l/m% correction in traditional non- relativistic expansions 141 ] . 

Extensive tests of the SR expansion have been performed within the context of the Fermi 
gas. In particular, comparisons with the fully relativistic Fermi gas (RFG), where the exact 
result is well known, have shown the reliability of the expansion |l[ • Similar SR expansions 
have also been performed for meson-exchange currents jl5j ]. nucleon- A currents [if]] and 
more recently for the charge changing (CC) weak current driving QE (vi,l~) reactions 

Part of the interest in reactions involving neutrinos instead of electrons lies in their 



implications for ongoing and planned neutrino oscillation experiments 
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201, 121|, [22|, 



26] . Since these naturally involve heavy nuclei as targets, reliable nuclear models 
of the reaction A(ui,l~) play an essential role in the interpretation of the data. Because of 
the close relationships that exist amongst all semi-leptonic electroweak processes, accurate 
descriptions of A(e, e') data appear as a requirement that sets strong restrictions for nuclear 
modeling of neutrino reactions. Based on what is known from previous approaches to this 
problem, and on the results presented in this paper, two ingredients arise as essential if the 
nuclear reaction modeling is to be successful in describing the electroweak cross sections in 
the kinematical regime of interest, namely, relativity and final-state interactions (FSI). 

First is relativity, since one is dealing with momentum transfers in the intermediate- 
to-high energy regime, typically of the order of 1 GeV or higher, for which tradi- 
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tional non-relativistic expansions in powers of 1/mjv are not applicable. Different fully 
relativistic models (based on Dirac equations and/or relativistic many-body theories) 
lave been developed in recent years aiming to describe electron and neutrino scattering 
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411 ]. While relativistic approaches 



based on the Dirac equation appear as the most direct way to deal with the problem, the 
importance of the different relativistic ingredients is more easily explored by detailed com- 
parison with the extensively employed non-relativistic approaches, and in particular to the 
SR one that is the focus of this work. This comparison will also help in identifying which 
relativistic ingredients (current operators, initial- or final-state wave functions) are the main 
ones responsible for the difference between relativistic and non-relativistic results. Upon 
incorporating these ingredients in SR approaches, one can extend their applicability to re- 
gions well into the relativistic domain, and thereby hope to produce reliable results at high 
energies. 

Secondly, the model of FSI must account for many-body effects that are known to be es- 
sential in describin g re asonably well (e, e') cross sections within the context of non-relativistic 



approaches 
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481 ] . For instance, medium modifications of the one-particle 



one-hole (lp-lh) final-state self-energy via two-particle two-hole (2p-2h) intermediate states 
in the continuum appear to be important in non-relativistic many-body treatments of a va- 



riety of processes involving low-to- intermediate momentum transfers |49l. 1501. \5l\. The main 
effect is a significant shift of strength to high energies in the QE cross section. On the other 
hand, relativistic impulse approximation modeling of inclusive reactions 4jJ, |52j shows an 
important effect of the FSI where strength is shifted to the high momentum region. 

The importance of FSI for intermediate-to-high momentum transfers has also been inves- 
tigated in connection with the scaling prop erties of QE cross section. The recent analysis 
of (e, e') world data in the QE region [53j, |5J, [55( has permitted the extraction of a uni- 



versal scaling function from the longitudinal response function data (here ip' is the 

dimensionless scaling variable defined in eq. (1291) below). The experimental scaling function 
presents an asymmetric shape, in contrast to the symmetric behavior predicted by most 
independent-particle models, which essentially give results that are similar to the RFG, 
with the exception of the tails observed for ip' < —1 and ip' > 1, where the RFG result is 
zero by construction. The behavior exhibited by the experimental scaling function reflects 
that of the longitudinal response function including the FSI for intermediate values of the 
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momentum transfer 



45| . i.e. the single-particle strength at the QE peak is reduced and a 



large tail is observed for high values of the energy transfer uj. Note that the difference with 
respect to previous approaches is that the experimental data now correspond to relatively 
high values of the momentum transfer, mainly in the range from q = 500 MeV/c to 1 GeV/c, 
for which scaling occurs, and where most of the existing non-relativistic models of FSI are 
no longer ap plic able. 



Recently 



41 



521 ] good agreement between the calculated scaling function and the data 



has been found in a fully relativistic approach based on the relativistic mean-field (RMF) 
model, within the impulse approximation (IA). The bound and outgoing nucleon states are 
described by the same self-consistent Dirac-Hartree potential. Here the use of the same 
relativistic potential in the initial and final states appears to be essential, since, as shown 
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52], calculations using the real parts of complex relativistic optical potentials do not 



produce the asymmetric behavior seen in the data. Apparently the main factors responsible 
for the asymmetry are not only relativistic kinematics, but also the particular dynamics 
contained in the Dirac equation when strong scalar and vector potentials are used for the 
continuum nucleon states (i.e., for the FSI). In fact, within the IA the significant shift of 
strength to positive ip' values only occurs under the presence of these strong potentials (41 ]. 

In this paper we investigate whether it is possible to describe the asymmetric behavior 
of the QE scaling function within the semi-relativistic approach for intermediate-to-high 
momentum transfers. As FSI are essential in order to perform meaningful comparisons with 
QE scaling data, we explore various ways of including these ingredients in the SR model 
at relativistic energies. We consider two different models of FSI: the Smith- Wambach 2p- 
2h damping (SWD) model and the Dirac-equation-based (DEB) potential plus the Darwin 
factor (DEB+D) 
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The SWD model is an extrapolation to relativistic kinematics of the nuclear re- interaction 



model introduced in 



43 



49], which gives one a straightforward way to incorporate the 



effects of 2p-2h excitations in the nuclear response function within a non-relativistic context. 
Speculations on how to extend this model to QE neutrino scattering for relativistic energies 



have been presented in 
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5911 . within the context of the RFC 



The DEB+D method attempts to translate the success of the RMF model in describing 



the phenomenological scaling function from QE (e, e') data 



41 



52j to the SR approach. 



We perform calculations in a SR continuum shell model where the final wave functions are 
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obtained by solving the Schrodinger equation with the DEB potential, and multiplying by 
the Darwin factor. The DEB potential is obtained by a reduction of the Dirac equation to a 
Schrodinger equation with a local potential, which implies that one must multiply the wave 



function by a non-locality Darwin factor 30 



56 



57|, 



601 ] . This essentially amounts to using 



upper components from the solutions of the Dirac equation in computing the final nucleon 
wave functions. 

Both the DEB potential and Darwin factor are functions of the local (energy indepen- 
dent) vector (V) and scalar (S) components of the relativistic Hartree potential. The non- 
relativistic reduction implied in the derivation of DEB potential and Darwin term introduces 
a linear dependency on the energy of the particle. The results obtained in the SR model de- 
scribing FSI in this way compare well with those obtained using the fully relativistic RMF 
model, and thus also reproduce the successful comparison with the experimental scaling 
function data. This allows us to conclude that it is the treatment of FSI in the RMF that 
gives rise to the large asymmetric tail in the superscaling function. 

In this paper we also present an application of the SR-DEB+D approach to the super- 
scaling analysis (SuSA) of neutrino CC QE cross section. This method has been proposed 
in [6j| as an efficient way of predicting neutrino cross sections from the (e, e') data, by ex- 
ploiting the scaling properties of the latter. The method is based on the hypothesis that 
the neutrino cross sections universally scale in the same way as do the electron scattering 



cross sections 
models 
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'or intermediate-to-high momentum transfers, as is seen for several types of 



52]. An application of the SuSA approach to compute integrated neutrino 



cross sections has been reported recently in 62|. In this work we investigate the validity of 
the superscaling approach within the SR-DEB+D model. Given the success of our model in 
describing the experimental scaling data, this check will help further to lay the foundations 
of the SuSA approach introduced originally in [61 1. 

The structure of the work is the following. In sect. II we briefly outline the SR model 
and the different treatments of the FSI. We present results in sect. Ill and our conclusions 
in sect. IV. 
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II. THE SEMI-RELATIVISTIC (SR) MODEL 

In this section we summarize the basic formalism for electron and neutrino reactions 
within the SR approach, including the different treatments of the FSI. We refer the reader 



to [1] and 



17J for specific details. 



A. Electron and neutrino cross sections 

We focus specifically on the reactions (e, e') and (v h l~) induced by electrons and neutri- 
nos, respectively, where l~ is a lepton with mass mi (typically a muon). The four-momentum 
of the incident lepton is k^ = (e, k) while k!^ = (e', k') is the four-momentum of the final 
lepton. The four- momentum transfer is denoted = k^ — k /fl = (u, q). The results below 
are referred to a coordinate system where the z-axis points along q and the x-axis along 
k - (k • q)q/g 2 . 

The inclusive cross section for these reactions can be written in general as 

where cq is the usual Mott cross sections for (e, e') reactions |4| , while it becomes an analo- 
gous factor in the case of <„, n re ac tions ( S ee QQ fol ^ ex presslonS ). T h e relevant 
observable is the nuclear structure function JF 2 which can be written as 



T 2 = v L R L + v T R T (2) 

for (e, e') and 

T 2 = VccRcc + 2VclRcl + VllRll + V t Rt + 2V T 'Rt' (3) 

for (v^l^) reactions. The lepton kinematical factors Vk and Vk are defined in Q] and 61], 



respectively. Note that, in the limit of lepton mass mi = 0, one has the identities vl = Vcc 
and vt = Vt- Finally, the nuclear response functions are given by 

Rl = Rcc = W 00 (4) 

Rcl = -\{W^ + W™) (5) 

Rll = W 33 (6) 

R T = W 11 + W 22 (7) 

Rt' = -\ {W 12 - W 21 ) , (8) 
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where the inclusive hadronic tensor is 

= J26(E f -E t - oo)(f\nQW(f\r(Q)\i). (9) 

fi 

Here J^(Q) is the nuclear current operator relevant for the reaction, i.e., the electromagnetic 
current in the case of (e, e') or the CC weak current in the case of (z/;, as specified below. 
Note that although in eq. (T4|) Rl and Rcc are formally equal, in practice they are not the 
same quantity, since different current operators and nuclear matrix elements are involved in 
their definitions, Rl referring to electron scattering and Rcc to neutrino reactions. 

The current operators used in this work are SR expansions of the fully relativistic ones 
in powers of r\ = p/mjv, where p is the momentum of the initial (bound) nucleon, which is 
typically small (r/ < 1/4), and therefore an expansion at most to first order (0(rj)) should 
be adequate for the inclusive reactions considered here in the region of the QE peak. The 



expansion was performed in [l| and [17J (see 4| for a review on the general SR expansion 
procedure). Here we give just the final expressions used in our calculations. The CC weak 
current is written in momentum space as J^ 1 = Jy — J A , in terms of the vector and axial 
current terms. For the vector current one has 

J°v = £o + *£o (KXr7)-<r (10) 

Jv = Zir^ + itiaxK (11) 

and for the axial current 

J° A = CS«.<r + CSV-o- (12) 

J' A = ('zK-v + CZ^-o- (13) 

Ji = C( ct\ (14) 

In eqs. ( flWl4j) we have introduced the dimensionless momentum transfer vector k = q/2mAr, 
while we use the notation J to denote the component of the vector J perpendicular to the 
momentum transfer q. 

Finally, the nucleon form factors and relativistic correction factors are included in the 
coefficients £j, Q and Q' (the corresponding relativistic versions of these quantities for 



8 



the electroweak neutral current were introduced in the appendix of pj), defined by: 
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(15) 
(16) 
(17) 

(18) 
(19) 



Here use has been made of the dimensionless variables k = q/2rriN, A = uj/2mN, and 
t = k 2 — A 2 . The nucleon form factors appearing in eqs. (fl5l - [T9l) are the isovector magnetic 
G V M = G P M — G\j, isovector electric G\ = G P E — G E and the axial form factor Ga- In Q 
and we have introduced the following combination of axial-vector and pseudoscalar form 
factors 

G' A = G A ~ tG p , (20) 

where Gp is the pseudoscalar nucleon form factor. 

The expression for the SR expansion of the electromagnetic current is very similar to 
that for the vector current in eqs. (11011111) . The only difference is that the proper proton or 
neutron form factors, instead of the isovector ones, should be inserted in the corresponding 
£i and ^ coefficients. 



B. The nuclear and FSI models 



In this work we describe the nuclear structure with an uncorrelated shell model. The 
many-body initial nuclear state \i) is a Slater determinant. It is constructed with a set of 
single-particle solutions of a non-relativistic mean-field potential. The final states |/) = 
are particle-hole excitations of the nuclear core. The holes \h), with non-relativistic 
energies e^, are solutions of the Schrodinger equation with a Woods-Saxon potential. The 
particles, \p), with asymptotic kinetic energies e p = + u, are continuum solutions of the 
Schrodinger equation with positive energy. In the SR approach the relativistic kinematics 
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is implemented by the following substitution of the eigenvalue for positive energies: 



1 



2m 



N 



(21) 



For intermediate energies is small compared with u and the above substitution can be 
replaced to a good approximation by A — > A(l + A) with small error [63]. 

We begin at the basic starting point, denoted below as semi-relativistic Woods-Saxon 
(SR-WS), where the continuum states are described by solutions of the Schrodinger equation 
with the same Woods-Saxon potential as the bound states. Building on this, two further 
different approaches to FSI have been considered. 

In a first treatment, we include the FSI effects coming from 2p-2h intermediate states 
by following the Smith- Wambach damping model (SWD). Within this model the total re- 



sponse function (q, to) is computed from the bare SR-WS one, R 



SR-WS 
K > 



or convolution integral 

POO 

R S k WD ^)= / dER s K R - ws (q,E) 
Jo 



m N 
M*(q) 



m N 
M*(q) 



E,-u 



folding 



(22) 



where M*(q) is the nucleon effective mass taken from 64J and, following |49j, the function 
p(E, uj) is given by 

1 Vd.iMO. 

(23) 



- 1 r H/ 2 
P( ' U) ~ 7T [E - u - A(cu)f + [r(a;)/2] 2 ' 



EM = AM 



4 r M 



(24) 



being the complex self-energy of the final states. The imaginary self-energy is computed 
through the following average for the particles and holes 



i r 

T(u) = - de [7 P (e) + j h (e - u)] . 
. — . ^ Jo 



We use the parametrization of 



65) 



7 P (e) = lh{e) = 7(e) = 2a 



e 2 + e 



(25) 



(26) 



with a = 10.75 MeV, e = 18 MeV and e x = 110 MeV. 

The second approach to FSI, and main focus in this work, is the Dirac-equation-based plus 



Darwin term (DEB+D). The DEB potential, Udeb-, is obtained from t 
written as a second-order equation for the upper component ip up ( r ) [30, 



re Dirac equation, 
""fil, 57)]. This wave 



function is then written in the form 



Vv(r)=^(r,£)0(r), 



(27) 
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and the Darwin factor K(r, E) is chosen in such a way that the function <p(r) satisfies a 
Schrodinger-like equation 

1 



2m 



N 



-V 2 + U DEB (r,E) 



E 2 -m 2 N 



2m 



N 



(28) 



Here E is the relativistic energy of the final nucleon. Note that both the DEB potential, 
Udeb{t, E), and the Darwin factor, K(r,E), show an explicit energy dependence coming 
from the non-relativistic reduction of the Dirac equation. The complete expressions in terms 
of the vector and scalar parts of the relativistic potential are given, for instance, in 30_|, l60 |. 
In the results section we show the energy dependence of these quantities in more detail. 



C. Superscaling approach 

In this paper we focus on the properties of the scaling functions /zXVO and /t^')- 
these are defined as the electromagnetic responses R^q^uj) and Rt(q,uj) divided by the 
appropriate single-nucleon functions, Gk {K = L,T), weighted by the nucleon number 



involved in the process (see 17 



61| for details). 

The scaling functions depend generally on q and u and are different for different nuclei. 
Scaling of first kind is said to hold when they are found to depend, not independently on 
q and u>, but only on a specific function of these two, namely the scaling variable ijj'(q,uj) 
defined as 



* = r F 1/2 - — — 2 , (29) 



;i + A')r' + /c v / H 1 + T ' 
where A' = (u — E s )/2m Nl k = q/2m N) t' = k 2 — A' 2 , and £p = a/1 + (hp/uiN) 2 — 1. 
Experimentally, the Fermi momentum kp and the energy shift E s are empirical parameters 
determined through fits to QE electron scattering data. The unshifted scaling variable ip is 
defined by the same formula, but with E s = 0. 

If fi and fx are independent of kp as well, one says that second- kind scaling in fulfilled. 
When both kinds of scaling occur, the responses are said to "superscale" . 

In the case of neutrino reactions, scaling functions can be obtained in a similar way, 



dividing the weak responses by the corresponding single-nucleon contribution, see 



17 



41 



611 ]. The approach called SuSA assumes that the scaling functions entering both (e, e') and 



(yi, I ) reactions are the same: in this case one can reconstruct the weak cross sections from 
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the electromagnetic ones by multiplying the scaling function by the weak single-nucleon 
contribution. 



III. RESULTS 

In this section we present results for the nucleus 12 C. The ground state is described here 
in the extreme shell model, namely the lsi/2 and l]?3/2 shells are fully occupied for protons 
and neutrons. The Woods-Saxon potential parameters used to describe the (non-relativistic) 
bound energies and wave functions are given in [l^j]. As stated in the previous section, for 
the final states in the SR model we use different forms of the potential depending on the 
description of the FSI. In the SR-WS the same Woods-Saxon potential is used for initial and 
final states. In the SR-DEB the DEB potential and the Darwin term are used for the final 
states. A general multipole expansion of the current operators is performed to compute the 
nuclear response functions and scaling functions. We refer the reader to H] for details 
on this aspect of the calculation. 

We first investigate the properties of the electromagnetic (e, e') L and T scaling functions, 
leaving the discussion of the neutrino CC cross section for the end of this section. In 
particular, we are interested in the study of the properties of scaling of the first kind. To 
this end we compute the function /^(^O {K = L,T) for fixed values of the momentum 
transfer q. One has scaling of the first kind when no dependence on q is observed. 



A. Relativistic effects 



We start with a brief discussion of the size of the relativistic effects embodied in the SR 
model. In Fig. 1 we show the longitudinal and transverse scaling functions of 12 C computed 
in the shell model for three values of the momentum transfer, q = 0.5, 0.7 and 1 GeV/c. 
We show results for three models of the reaction, all of them including the same Woods- 
Saxon potential in initial and final states. The dashed lines correspond to a traditional 
non-relativistic approach using non-relativistic currents and non-relativistic kinematics. The 
dotted lines have been obtained using relativistic kinematics, but still the non-relativistic 
current operators. Finally, the solid lines correspond to the SR approach with relativistic 
corrections in the current as well, i.e., to the SR-WS model. 
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FIG. 1: Scaling functions of C in the continuum shell model with a Woods-Saxon potential. Solid 



lines: SR model. Dashed lines: Traditional non-relativistic results. Dotted lines: Non-re 
current operators using relativistic kinematics. The experimental data are taken from 



ativistic 



As one can see, the width of the scaling functions is significantly reduced in the calcu- 
lations that include relativistic kinematics, while the inclusion of relativistic effects in the 
current operators leads to an increase of fi and a reduction of fx- The importance of rela- 
tivistic effects clearly increases with the momentum transfer. In the figure we also plot the 
experimental data for the averaged /lC0')> taken from 55]. It is important to point out 
that data refer only to the analysis of the longitudinal scaling function, i.e., fi; however, we 
have chosen to present such data also in comparison with theoretical results for fx- This 
makes the comparison between fi and fr more straightforward and, in fact, it allows us to 
check the degree to which scaling of the zeroth kind, i.e., fi = fr, is fulfilled for the different 
models presented in this work. 

Incidentally, note that the non-relativistic current with relativistic kinematics (dotted 
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FIG. 2: In each panel the 12 C scaling function is displayed for three values of q = 0.5 (solid 
lines), 0.7 (dashed) and 1 GeV/c (dotted). The initial and final states are described with the 
same Woods-Saxon potential. Upper panels: Traditional non-relativistic model. Middle panels: 
Non-relativistic model using relativistic kinematics. Lower panels: SR model. 



lines) gives the best description of experimental data for q = 1 GeV/c, where a non- 
relativistic treatment of the current operator should hardly be adequate. This fact is a 
consequence of the relativistic factor n/y/r in the coefficient £ in eq. (|T5|) which produces 
a reduction of Jl- 

Although from the full set of results seen in Fig. 1 it clearly appears that none of the 
three approaches provides a satisfactory description of the experimental data, one can still 
study the scaling properties of these models. One knows from previous superscaling studies 
that scaling of the first kind is approximately valid for the experimental longitudinal scaling 
function, f^. This places important restrictions on theoretical models, namely, that they 
should be (reasonably) compatible with scaling of the first kind. Otherwise they would be 
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unable to describe the (e, e') cross section over a wide range of momenta. This is illustrated 
in Fig. 2 for the different approaches in the shell model. There we show the same results as in 
Fig. 1, but presented in a different way, namely, this time plotting together the results from 
each of the models for the three values of q simultaneously. From Fig. 2 one can appreciate 
that the SR approach is the only model that shows scaling of the 1st kind. Therefore, we 
conclude that the relativistic corrections that are consistently accounted for in the SR, in 
both the current operator and the kinematics, are important for recovering scaling of 1st 
kind, although the model still lacks dynamical ingredients from relativistic FSI needed to 
reproduce the shape of the experimental data, as discussed below. 

Before entering into a detailed discussion of FSI, in Fig. 3 we display a test of the SR 
expansion, by comparing with a fully relativistic result. Similar tests have been carried out 
in the context of the relativistic Fermi gas in l[ However, at the level of the SR-WS it 
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is not possible to make meaningful comparisons with fully relativistic results, because of the 
different description of the dynamics in the final states. To make a consistent comparison 
one should go a step back and consider the plane- wave impulse approximation (PWIA), 
where the final potential is set to zero. This is done in Fig. 3, where we show the PWIA 
results for the L and T scaling functions in semi-relativistic and relativistic models. From 
these results one can see that the description of the initial nuclear state and of the current 
operator are quite similar in both models. Note that in the relativistic calculation the CC2 
current has been employed (see also later). The bound wave functions and energies are very 
close in the relativistic and non-relativistic approaches. The small differences seen in Fig. 3 
are linked to off-shell effects in the current matrix elements, and are of the same order as can 
be found between relativistic results using different prescriptions for the current operator 
41]. 

B. FSI in the damping model 

There are two main effects embodied in the Smith- Wambach damping model used here: a 
shift and a redistribution of strength, as illustrated in Fig. 4, where we show the longitudinal 
response function for five values of the momentum transfer. These effects are a consequence 
of the folding integral in eq. ( |22i) . which basically produces the redistribution, and of the 
appearance of the effective mass M*(q), which is responsible for the shift. However, when 
the energy u is large, the parametrization of eqs. (I25II26I) results in a small nucleon width: 
hence the Lorentzian function p(E, uj) becomes close to a Dirac delta function and only the 
effective mass effect remains. This is why in Fig. 4, for the largest g-value, q = 1.5 GeV/c, 
the FSI just produce a shift of Rl- Thus, although the damping model can produce the 
needed effect for low-to-intermediate momentum transfers (i.e. redistribution of the strength 
as required by the data), for higher momentum transfer this method is basically equivalent 
to introducing an almost constant effective mass (of around M* /mjy — 0.8 — 0.9). 

The scaling properties of Rl under the SWD model of FSI may be seen by examining the 
results in Fig. 5. There one observes that the scaling function increases with q, and hence 
that scaling of the first kind is not respected by the model. Moreover it is almost symmetric, 
whereas the experimental data, also shown in the figure, display an asymmetry, indicating 
stronger re-distribution of strength for positive values of ip' . The apparent incapability of the 
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FIG. 4: Effect of FSI on the longitudinal response function. Dashed lines: results in the SR model 
with a Woods-Saxon potential. Solid: including FSI in the Smith- Wambach damping model. 
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FIG. 5: Longitudinal scaling function in the SR model with Smith- Wambach damping. Several 
values of the momentum transfer q = 0.5 (solid lines), 0.7 (dashed), 1.0 (short-dashed), 1.3 (dotted) 
and 1.5 GeV/c (dot-dashed) are plotted. No energy shift has been applied to the theoretical 



calculation. Experimental data from 



55] 



SWD to describe this feature of the data should not be considered a failure of the model, but 
simply an inadequate extrapolation to high energy of a model that was originally proposed 
to describe the damping of the continuum nuclear response in the region of small energy 
transfer. 



C. The DEB approach to FSI 

The second approach to FSI in this work amounts to using the DEB potential plus 
Darwin term as described above. The interest in this approach for intermediate-to-high 



energies originates from the results found in 



41 



521 ] where the relativistic mean-field (RMF) 



model has been found to describe reasonably well the scaling function data. In particular, 
the asymmetry presented by the data appears to be obtained only in IA-based models when 
one uses the same form of the (real) S and V relativistic potential in both initial and final 
states. One must at this point remember that the present focus is on inclusive reactions: one 
should not be tempted to apply the optical potential fit to elastic data to inclusive reactions 



such as (e, e') or (z/;, I ), where different channels are open 66(. It is not the purpose of this 



paper to find rigorous theoretical arguments to justify the use of the Hartree potential as 
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FIG. 6: Effect in the DEB+D model of FSI on the L and T responses. Dashed lines: SR-WS 
model. Solid: SR-DEB+D. Dotted: SR-PWIA. 



the best choice in relativistic calculations. Developments along these lines will be discussed 
elsewhere 671 ] . We just mention here that the asymmetry found in the scaling function 
computed in the RMF is a consequence of the large scalar and vector parts of the potential 
used in the final state, which are both much smaller in the case of the real parts of typical 
optical potentials or in other similar approaches. 
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The goal of this section is to show that when conveniently implemented in the SR ap- 
proach, the same relativistic Hartree potential yields results consistent with the ones ob- 
tained within the fully relativistic mean-field RMF model. The electromagnetic response 
functions with the DEB+D choice for FSI are shown in Fig. 6 for fixed momentum transfer 
in the range q = 0.5 to 1.5 GeV/c. In the figure we show the SR results using the Woods- 
Saxon potential in the final state (dashed lines) and using the DEB+D model of FSI (solid 
lines). For comparison in the same figure we also show the SR-PWIA results. As one can 
see, the effect of the DEB+D model is precisely a shift and re-distribution of the strength, 
typically what one expects to occur due to FSI mechanisms. For q = 0.5 GeV/c the effect is 
similar to the SWD results of Fig. 4. However, for higher values of q the effect of the FSI is 
maintained in the T response and even increases in the case of the L response, contrary to 
what happens in the SWD model. In this way, for q = 1.5 GeV/c the total strength is not 
only re-distributed, but amplified in the longitudinal response function. The behavior of Rl 
for high q in the high-tu tail is just a consequence of the SR expansion used for the current 
operator. In fact the time component of the vector current in eq. ( fTUl) is proportional to the 
£o coefficient which in turn is proportional to the relativistic factor k/\/t in eq. (115p . The 
fact that this factor becomes infinite on the light cone u = q is just a consequence of the 
off-shell properties of the SR expansion, since it was performed for on-shell nucleons and 
therefore its applicability should be most appropriate for the QE-peak region. This is not a 
fault of the current, since similar anomalous off-shell effects can also be seen in prescriptions 



of relativistic current 



411 ] . Therefore, one should be careful when extrapolating these oper- 



ators to highly off-shell conditions, and the results of Fig. 6 should be taken with caution 
in the region close to u = q, in particular when q is high and non-QE effects dominate the 
total cross section, as there probably Rl cannot be extracted experimentally. In fact the 
experimental data for the scaling function shown in this paper were taken from world data 
in the range q = 0.5 to 1 GeV/c. 

The interest in showing results for values of q as high as 1.5 GeV/c is related to the study 
of the properties of scaling of the first kind performed next. In Fig. 7 we show the scaling 
functions in the DEB+D model. They have been extracted from the results of Fig. 6, and 
are displayed as functions of the unshifted (ip) and shifted {if}') scaling variables. When 
nonzero, the energy shift is 20 MeV in all cases. The experimental data for fi are shown 
for comparison. One can see that the DEB+D results follow the trend of the experimental 
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FIG. 7: First-kind scaling in the DEB+D model. In each panel results are shown for q = 0.5, 
0.7, 1.0, 1.3 and 1.5 GeV/c. First row of panels: scaling functions versus ip, namely, when E s = 0. 
Second row: the same versus ip' with constant energy shift E s = 20 MeV. Experimental data from 
55]. 

data, clearly showing the same asymmetric shape. While the strength and width of the 
results of Fig. 7 are basically the same as the data, the theoretical results are shifted to the 
right with respect to the data. Though small, the shift increases with q. Due to the shift, 
one would conclude that scaling of the first kind is not perfect in this model. However, note 
that the data do not contain values of q as large as 1.5 GeV/c and that the error bars are 
a consequence of a small dispersion, since scaling of the first kind is not experimentally as 
good as scaling of the second kind (i.e., independence of the nucleon number A). Therefore 
a small scaling violation as that seen in Fig. 7 may not actually be inconsistent with the 
data. Another interesting conclusion from these results is that to a high degree the DEB+D 
model respects scaling of zeroth kind, i.e., fi = fr, for all q- values. 

Another study of scaling of the 1st kind within the DEB+D model has been performed 
with results at a different kinematical setting shown in Fig. 8. This time we show the 
scaling functions computed for several incident electron energies and fixed electron scattering 
angle 9 e = 45°. Here the momentum transfer is not constant and changes with ip. As a 
consequence, the scaling works better and the curves are closer to each other. 

At this point a test of the present SR model is possible by comparing with a fully rela- 
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FIG. 8: First-kind scaling in the DEB+D model, 
electron energies e e = 1, 1.5 and 2 GeV. The elect 
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In each panel the results are shown for incident 
on scattering angle is 9 e = 45°, and E s = 0. 



tivistic model. Since the final states jn the DEB+D approach are computed essentially as 
in the RMF approach presented in 
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521 ]. both models should give similar results if the 



relevant relativistic corrections are included properly in the SR approach. A first compari- 
son is performed in Fig. 9 for three values of the momentum transfer. The RMF results are 
computed with the CC2 current operator. The SR model using the Woods-Saxon potential 
in the final state is also shown for comparison. It is apparent from the figure that the be- 
havior of the SR-DEB+D results is similar to that in the fully relativistic approach. The 
differences seen between the two models can partially be attributed to the peculiar off-shell 
properties of the CC2 and SR currents, as well as to the four-component structure of the 
wave functions involved in the RMF model compared with the SR-DEB+D. In fact differ- 
ences of the same order can also be found in the RMF approach between results obtained 
using other current prescriptions. An example is shown in Fig. 10, where we compare our 
results for q — 1 GeV/c with the RMF results employing two current prescriptions. For this 
value of the momentum transfer, the DEB+D model gives similar results to those for the 
RMF model with the CC3 current operator, although not exactly the same. From the results 
of Figs. 9 and 10 we conclude that our model contains the relevant relativistic corrections 
in this kinematical region. Our current operator could then be considered at this level as 
another prescription for off-shell extrapolation of the relativistic nucleon current. At this 
point it is also interesting to note that the RMF model with the CC2 operator leads to a 
clear violation of zeroth-kind scaling (likewise for CC1, but not for CC3), contrary to what 
happens with SR-DEB+D. 
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FIG. 9: Electromagnetic scaling functions for three values of the momentum transfer. Dotted lines: 
SR model with a Woods-Saxon potential in the final state. Solid: SR model using the DEB+D 
approach for FSI. Dashed: RMF with the CC2 current operator. Here E s = 0. 



Motivated by the results of Fig. 7, where scaling of first kind is fulfilled except for a small 
shift, and in order to explore further the scaling properties of our model, we have performed 
the calculation of the scaling variable ip' using a g-dependent energy shift E s (q). The value 
of E s (q) is fitted to get the maximum of J'l at ip' = 0. The results of this fit are shown in 
Fig. 11. One sees that the collapse of the curves is improved by this procedure. Although 
the scaling is not perfect, it is remarkable that for a wide range of q values (q = 0.5 to 1.5 
GeV/c) our results give essentially the same scaling function. The scaling is better for the 
transverse function f T , while the width of fi slightly increases with q. In Fig. 11 we also show 
a comparison with the experimental data. The agreement between theory and experiment is 
excellent. We recall that data refer only to the analysis of the longitudinal response, hence 
caution should be exercised when comparing the calculated transverse scaling function 
and data. Finally, we also present in Fig. 11 the shift energy E s (q) as a function of q, 
showing that the dependence of the fit is clearly linear. This linear behavior of E s (q) is 
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FIG. 10: Electromagnetic scaling functions for q = 1 GeV/c. Solid lines: SR model within the 
DEB+D approach for FSI. Dashed: RMF with the CC2 current operator. Dotted: RMF with the 
CC3 current operator. 



connected to the energy dependence of the DEB potential, which is also linear [30, |60j • We 
illustrate this dependence in Fig. 12, where we plot the central and spin-orbit parts of the 
DEB potential for a range of final nucleon kinetic energies from 100 to 1500 MeV. More 



details on the properties of the DEB potential can be found in [30j. From Fig. 12 it is 
evident that the central part of the DEB potential increases linearly with the energy and 
that it is very repulsive for high energy. The spin-orbit part, on the other hand, is rather 
small and unimportant for these energies. In the same figure we also show the Darwin factor 
K(r, E) for the same nucleon energies. 

The highly repulsive behavior of the DEB potential is the main reason for the strength re- 
distribution seen in our results. A repulsive potential favors the emission of nucleons having 
high enough available energy. Accordingly, more strength is placed in the high energy tail of 
the cross section or, equivalently, at positive values of ip' in the scaling function. One should 
mention here that the DEB potential alone is not enough to produce reasonable results for 
the scaling function. The Darwin factor, K(r,E), shown in the last panel of Fig. 12, must 
also be applied to the wave function. We illustrate this point in Fig. 13, where we compare 
the SR-WS model (namely with Woods-Saxon potential in the initial and final states) with 
the SR-DEB, but with the Darwin term set to one. One can see that the DEB potential 
alone produces the mentioned shift and widening of the scaling function to higher energies, 
due to the repulsive character of the potential, although the strength is still too high when 
compared with the data. 

The effect of the Darwin term can be appreciated in Fig. 14, where we plot the scaling 
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FIG. 11: First-kind scaling in the SR-DEB+D model. Results are shown for q = 0.5, 0.7, 1.0, 
1.3 and 1.5 GeV/c. First row: Scaling functions versus ip' using a ^-dependent energy shift E s (q). 
Second row: the same compared with the experimental data. Last panel: the linear behavior of 
E s (q). The fitted values of E s are 20, 32, 50, 68 and 80 MeV. Experimental data from 55]. 



functions in the SR-DEB approach with and without that factor (or equivalently by setting 
K(r,E) = 1 in the calculation). Since K(r) < 1 inside the nucleus, its effect is a reduction 
of the scaling function. This reduction is precisely the one needed to reach the experimental 
data, also shown in Fig. 14. Thus the Darwin factor is essential in this approach to FSI. 

One should mention that the full DEB+D model is not equivalent to solving a Schrodinger 
equation with a Hermitian local potential. Actually, the wave function in eq. (|27"1) is a 
solution of such an equation with a non-local potential, since apart from the dependence 
on the energy E, dependence on the momentum or gradient operator also appears from 



the reduction of the Dirac equation 30 



56 



57 



60(. This means that the solution of the 
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FIG. 12: First row: DEB potential for final nucleon kinetic energies 100-1500 MeV. Last panel: 
the same for the Darwin factor. 

Schrodinger equation with DEB potential plus Darwin factor is not comparable to any 
other non-relativistic or relativized solutions obtained with local potentials. The non-locality 
contained in the DEB+D approach introduced here, also implicit in the Dirac equation and 
made explicit by means of the DEB+Darwin reduction, thus becomes essential in producing 
the kind of effects that usually appear in non-relativistic approaches of FSI related to the 
role of correlations and/or exchange terms in the optical potentials. In particular, one could 
argue that in the SWD model discussed above, the inclusion of 2p-2h intermediate states 
gives rise also to non-local effects in the final state. 

A better appreciation of the effect of the DEB potential without the Darwin factor can 
be seen in Fig. 15, where we compare the scaling functions in the SR-WS and SR-DEB 
potentials. Once again the Darwin factor is set to one. The g-dependent shift E s (q) has 
been applied to the DEB results. The widening of the distribution to higher energies, 
producing a longer tail, is clear. Note also that the DEB potential does not modify the 
results in the region ip' < 0, and that there is a reduction of the maximum in the region 
ip' ~ (the magnitude of the reduction increases with q) . 
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FIG. 13: Effect of the DEB potential without the Darwin term. Solid lines: SR model with the 
same Woods-Saxon potential in the initial and final states. Dashed: SR-DEB results with the 
Darwin term set to one. Experimental data from 551 ] . 



D. Applications to neutrino reactions 

We now exploit the SR-DEB+D model introduced in the last section to make predictions 
for QE CC neutrino reactions. Given the success of that model in describing the (e, e') 
superscaling data, the neutrino cross section results obtained with our model can be con- 
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FIG. 14: Effect of the Darwin term on the scaling function in the SR model. Solid: total DEB+D 
result. Dashed: results with the Darwin term set to one. Experimental data from 551 ] . 



sidered at least reasonable in a kinematical range from intermediate-to-high energies and 
momentum transfers. 

A selection of results for the 12 C(f M , differential QE cross section are presented in 
Figs. 16 and 17, for incident neutrino energies of 1 and 1.5 GeV, respectively. We show 
results for a range of scattering angles from 15° to 150°. All of the results in the figures 
have been obtained with the SR model with different ingredients for the FSI. The results 
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FIG. 15: Effect of the DEB potential without the Darwin term. Solid lines: SR model with the 
same Woods-Saxon potential in the initial and final states. Dashed: SR-DEB results with the 
Darwin term set to one and with the ^-dependent energy shift E s (q). 



within the DEB+D model of FSI are shown with solid lines, while the dotted lines are the 
SR-WS results using the same Woods-Saxon potential in initial and final states. We display 
the cross section as a function of the final muon energy e'. Therefore the tail produced by 
the FSI can now be seen to the left of the peak of the SR-WS results. Finally, with dashed 
lines we show the results obtained in a third calculation which closely follows the general 
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FIG. 16: Double-differential cross section for the CfV^j H~) reaction for incident neutrinos of 
1 GeV, computed with several versions of the SR model. Dotted: Wood-Saxon potential. Solid: 
DEB+D. Dashed: reconstructed from the (e, e') DEB+D model using the SuSA. 



procedure originally introduced in 



611 ] and referred to as SuSA approach. 



Just as when the SuSA approach was performed using data, here the theoretical (e, e') 
cross section is divided by the single-nucleon factor to obtain the scaling function, which 
in turn is multiplied by the relevant CC weak factor, to obtain the neutrino cross section. 
That is, this procedure is similar to the SuSA approach in [6l|] except that instead of using 
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FIG. 17: The same as Fig. 16, but for neutrinos of 1.5 GeV. 



the "experimental" superscaling function obtained from the analysis of (e, e') world data, 
in this case we make use of the scaling function derived from (e, e') results with the SR- 
DEB+D model. Notice that the reconstruction of the (z/^/i - ) cross section from the (e, e') 
data (SuSA) or (e, e') calculations (we denote it as superscaling-based approach to avoid 
confusion with SuSA which implies the "experimental" scaling function) assumes that the 
scaling function for both reactions is similar for the kinematics of interest. Thus in Figs. 16 
and 17, we treat the solid lines of the DEB+D model as "exact" and check the accuracy 
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of the superscaling-based approach by comparison of the dashed lines with the solid ones. 
One can see that in general the superscaling-based results are very close to the DEB+D 
results. The differences between the two sets of results become smaller as the scattering angle 
and neutrino energy increase, since bigger values of the momentum transfer are involved. 
The biggest differences appear for very small angles, implying rather small values of the 
momentum transfer. In these cases the description of the giant resonance region needs a 
nuclear model containing collective degrees of freedom such as the RPA, not included in 
our approach. An example is shown for the case of 9 = 15° in Figs. 16 and 17. The 
presence of the narrow peaks appearing as potential resonances clearly violates the scaling 
of the cross section in that region. However, for higher values of the scattering angle the 
superscaling-based analysis begins to be applicable, and its predictions can be considered 
quite reasonable, with an error typically below 10%. 

IV. CONCLUSIONS 

In this work we have presented improvements over the semi-relativistic approach to elec- 
tron and neutrino reactions for intermediate-to-high energies and momentum transfers in 
the QE region. Going a step beyond the continuum shell model, we have explored two ap- 
proaches to describe the FSI: the Smith- Wambach damping model and the Dirac-equation- 
based potential. The former is applicable to low energies and has simply been extrapolated 
here to relativistic energies, while the latter is more appropriate for the kinematical regime 
considered in this paper and has been our main focus. 

Thus, in addition to using the semi-relativistic expansion of the nucleon current in powers 
of p/itln and relativistic kinematics, we have used the DEB form of the relativistic Hartree 
potential term to describe the continuum wave function of the ejected nucleon. Furthermore, 
we have included in the wave function the non-localities arising from the reduction of the 
Dirac equation by multiplying it by the corresponding Darwin factor. 

Firstly, we have focused on the analysis of the electromagnetic response functions. In 
particular we have investigated the properties of scaling of the first kind displayed by our 
results. The longitudinal and transverse scaling functions have been computed for a wide 
range of momentum transfers, from 0.5 to 1.5 GeV. Several aspects and details of the different 
theoretical ingredients embodied in our model have been analyzed. In particular, we have 
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examined the effects of the relativistic corrections and we have presented comparisons with 
fully relativistic results. 

We have found that our model approximately fulfills scaling of the first kind except for 
a small energy shift, E s (q), which turns out to be linear in q. This behavior has been con- 
nected to the repulsive character of the DEB potential, which also depends linearly on the 
nucleon energy. The Darwin factor has been found to be essential for the description of 
the experimental scaling function data. The Darwin factor is needed to correct for mathe- 
matically unavoidable non-localities arising in the differential equation describing the upper 
component of the relativistic nucleon wave function. 

Is is a remarkable result of this work that, except for the energy shift, our model gives 
essentially the same scaling function for a 1 GeV wide range of three-momenta q — spanning 
a kinematical region that extends from non-relativistic to relativistic conditions - - and 
that the behavior of the theoretical scaling function is essentially the same as that of the 
experimental data. The study of the theoretical energy shift performed here can be of help 
for future analyses of the scaling properties of experimental data. In particular, it suggests 
that the small scaling violations of the first kind found in the (e, e') data can be used to 
extract valuable information about the strength of the FSI. 

Finally, we have presented an application of the model to the inclusive CC neutrino 
reaction 12 C(z/ M , fi~) in the region of the QE peak. Results have been presented in a range 
of kinematics for several neutrino energies and muon scattering angles. In particular, we 
have used our model to investigate theoretically the kinematical range of validity of the 
superscaling approach in reconstructing the neutrino cross section from the electromagnetic 
scaling function. The validity of the superscaling-based approach has been verified with our 
model for kinematics involving values of the momentum transfer large enough such that the 
QE region is not contaminated by the presence of giant resonances. For incident neutrinos of 
1 GeV, this happens for angles typically bigger than ~15°. In these cases, within our model 
the superscaling-based approach is seen to be very successful in being able to reconstruct 
the neutrino cross section from the electromagnetic one. 
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